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SUMMARY 

The reliability with bounded distribution parameters (mean and standard deviation) was maxi- 
mized and the reliability-based cost was minimized for adaptive intraply hybrid fiber composites by 
using a probabilistic method. The probabilistic method accounts for all naturally occurring uncertain- 
ties including those in constituent material properties, fabrication variables, structure geometry, and 
control-related parameters. Probabilistic sensitivity factors were computed and used in the optimization 
procedures. For an actuated change in the angle of attack of an airfoil-like composite shell structure 
with an adaptive torque plate, the reliability was maximized to 0.9999 probability, with constraints on 
the mean and standard deviation of the actuation material volume ratio (percentage of actuation com- 
posite material in a ply) and the actuation strain coefficient. The reliability-based cost was minimized 
for an airfoil-like composite shell structure with an adaptive skin and a mean actuation material vol- 
ume ratio as the design parameter. At a 0.91 -mean actuation material volume ratio, the minimum cost 
was obtained. 


INTRODUCTION 

Aerospace structures are complex assemblages of structural components that operate under severe 
and often uncertain service environments. They require durability, high reliability, light weight, high 
performance, and affordable cost. Composite materials are potential candidates that meet these require- 
ments because they possess outstanding mechanical properties with excellent fatigue strength and cor- 
rosion resistance. Their mechanical properties are derived from a wide variety of variables such as 
constituent material properties and laminate characteristics (fiber and void volume ratios, ply orienta- 
tion, and ply thickness). These variables are known to be uncertain. 

A new challenge is to further enhance the structural performance by investigating other advanced 
concepts. Recent developments in smart structure concepts that use actuation materials such as piezo- 
electric ceramics have the potential to improve structural performance, durability, and reliability (refs. 1 
and 2). The control devices in smart structures consist of (1) a polarized material, (2) an electric field 
parallel to the direction of polarization, and (3) the expansion-contraction effects of the polarized 
material. When a control voltage is applied, the actuation material expands or contracts so that the 
structural behavior is altered by a desired amount and its reliability is changed. Present piezoelectric 
technology has been successfully applied to small-scale, low-stress structures. However, there are 
inevitable difficulties when current technology is applied to large-scale, high-stress composite 
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structures. These difficulties can be alleviated if the following concept is adopted: combine special 
piezoelectric fibers that have a fast actuation capability and relatively high-strength, high-modulus fiber 
to form an adaptive intraply hybrid composite. This composite can be readily integrated into a smart 
composite structure by using combinations of intraply and interply hybrid composites to ascertain that 
smart composite structures will operate in the design-specified range. 

The adaptation of the intraply hybrid composite concept (refs. 3 and 4) to smart composite struc- 
tures is depicted schematically in figure 1, which shows that the smart composite comprises (1) regular 
plies that consist of regular composite materials and (2) control plies that consist of regular strips of 
regular composite materials and strips of mixed regular and actuation materials. Actuators, containing 
actuation materials such as piezoelectric ceramics or piezoelectric fibers, are used to control the 
behavior of the composite structure by expanding or contracting the actuation strips to achieve the 
requisite design and operational goals. However, the strains induced by the actuator are affected by 
uncertainties in several factors that can only be quantified probabilistically. These factors include 
(1) inaccurate measurements made by sensors, (2) deviation from the intended electric field, (3) uncer- 
tain actuation strain-electric field strength relationship, (4) uncertain material properties for actuation 
materials, (5) uncertain electric field strength, and (6) improper location of sensor and/or control 
materials. Because of these factors, the use of control devices increases the uncertainty in already 
uncertain composite structural behavior. 

To account for various uncertainties and to satisfy design requirements, knockdown (safety) 
factors are used extensively. These knockdown factors result in a substantial weight increase but with- 
out a quantifiable measure of reliability. This paper describes an alternate method to determine struc- 
tural reliability. This method is embedded in the computer code IP ACS (Integrated Probabilistic 
Assessment of Composite Structures, ref. 2) for a comprehensive probabilistic assessment of composite 
structures. The schematic of IPACS is shown in figure 2. 

Because actuation materials are more expensive than regular composite materials, optimization 
techniques should be sought to achieve the balance between maximum reliability and minimum cost. 

In this paper, two types of optimization issues with reliability considerations are discussed: reliability 
maximization with parameter constraints and reliability-based cost optimization. 


SYMBOLS 


C F cost due to failure 

Cj manufacturing cost 

Cj(Pj) manufacturing cost for j th distribution parameter 

C 0 constant cost 

Cj total cost 

C a constant 

d 31 actuation strain coefficient 
E fll fiber modulus in longitudinal direction 
Ef 22 fiber modulus in transverse direction 
E m matrix elastic modulus 

G fl2 in-plane fiber shear modulus 

Gf 23 out-of-plane fiber shear modulus 
G m matrix shear modulus 

m ; mean of i* random variable failure 

nv, i lh reference mean 


2 



N number of distribution parameters 

P f probability of failure 

Pj i* distribution parameter 

Pj lower bound for Pj 

Pj upper bound for Pj 

Pj j* distribution parameter 

Pj mean value of actuation material volume ratio 

SF safety factor 

AT temperature difference 

u* most probable point 

V electric field strength 

Xj i* random variable 

Z structural performance 

a thermal expansion coefficient 

P reliability index 

A. electrically induced in-plane strain 

v fl2 in-plane fiber Poisson’s ratio 

Vj 23 out-of-plane fiber Poisson’s ratio 

v m matrix Poisson’s ratio 

o t standard deviation of i* random variable 

o Q reference standard deviation 

a 0 i* reference standard deviation 

O cumulative distribution function of normal random variable 

Subscripts: 

i, j indices 


IPACS COMPUTER CODE FOR RELIABILITY CALCULATION 

The IPACS computer code (ref. 5) has evolved from extensive research activities at the NASA 
Lewis Research Center to integrate probabilistic structural analysis methods (ref. 6) and computational 
composite mechanics (ref. 3). The composite micromechanics, macromechanics, and laminate theory 
(including interply and intraply hybrids) are embodied in ICAN (Integrated Composite ANalyzer, ref. 3). 
IPACS consists of two stand-alone computer modules: PICAN (Probabilistic Integrated Composite 
ANalyzer) and NESSUS (Numerical Evaluation of Stochastic Structures Under Stress). PICAN simu- 
lates probabilistic composite mechanics (ref. 7); NESSUS uses the information from PICAN to simu- 
late probabilistic structural responses (ref. 8). Direct coupling of these two modules makes it possible 
to simulate uncertainties in all inherent scales of the composite — from constituent materials to the 
composite structure, including its boundary and loading conditions as well as environmental effects. It 
is worth noting that the special reliability algorithm FPI (fast probability integrator, ref. 9) is used 
instead of the conventional Monte Carlo simulation (ref. 10) to achieve substantial computational 
efficiencies that are acceptable for practical applications. Therefore, a probabilistic composite structural 
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analysis which cannot be done traditionally becomes feasible, especially for composite structures that 
have a large number of uncertain variables. The results from the IP ACS analysis include a probability 
distribution function of the structural response, the reliability for a design criterion, and the probabilis- 
tic sensitivity factors of the primitive variables to the structural response and structural reliability. 

In IP ACS computational simulation for smart composite structures, electrically induced in-plane 
strain X is defined by 


X = d 31 V (1) 

where d 31 is the actuation strain coefficient, and V is the electric field strength. Because of the similar- 
ity between thermal strain and strain induced in the actuation material, the actuation strain is simulated 
using the thermal strain computed from a temperature field and the thermal expansion coefficient: 

X = d 31 V = a AT (2) 

where a is the thermal expansion coefficient, and AT is the temperature difference. 


PROBABILISTIC SENSITIVITY INFORMATION FOR OPTIMIZATION 

The commonly used sensitivity in a deterministic analysis is the performance sensitivity 3Z /3X i? 
which measures the change in the performance Z due to the change in a design parameter Xj. This 
concept is extended to a probabilistic analysis to define the probabilistic sensitivity that measures the 
change in reliability relative to the change in each random variable. The failure probability for a given 
performance is defined in the following equation (ref. 11): 


P f = <K-P) 


(3) 


where d> is the cumulative distribution function of a normally distributed random variable, and p is the 
reliability index. The probabilistic sensitivity factor (PSFj) for i* random variable is defined as 


PSF : 


3p _ u i 
dX: P 


(4) 


where u* is the most probable failure point of a limit state function in a unit normal probability space. 
The sensitivity of design parameters to structural reliability is another bit of useful information for 
controlling and adjusting design parameters from manufacturing to obtain the best benefit with mini- 
mum alteration. The sensitivity of the reliability to the mean of a normally distributed random variable 
Xj can be represented by the following equation (ref. 11) 

3L.-Hi < 5 ) 

3m j (J; 
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where nij and (jj are the mean and standard deviation of random variable Xj, respectively. Similarly, 
the sensitivity of the reliability to the standard deviation of a normally distributed random variable Xj 
can be computed from 

dfi _ _ SF, n* = _ uf (6) 

"5a" a, m 

With this information, optimization with reliability considerations can be achieved as demonstrated in 
the following sections. 


OPTIMIZATION WITH RELIABILITY CONSIDERATIONS 

Two optimization problems with reliability considerations are addressed and discussed in the 
following sections. 


Maximization of Structural Reliability with Bounded Distribution Parameters 

The constraints are applied to distribution parameters such as mean and standard deviation for each 
random variable as shown in the following equation: 

Pj £ Pj < P u j (7) 


where P; is the 1 th distribution parameter and P-, and P°j are the lower and upper bounds for Pj, respec- 
tively. From equation (5), it can be seen that the reliability can be increased by increasing the mean of 
the random variable with a negative probabilistic sensitivity factor and vice versa. In equation (6), it is 
found that 3(3/30, is always negative. Therefore, reliability can be increased only by scatter reduction 
regardless of the sign of the sensitivity factor. With these guidelines, reliability maximization can be 
easily accomplished. 


Minimization of Reliability-Based Cost 

A major design goal is to achieve the balance between maximum reliability and minimum cost. 
The criterion that simultaneously addresses both reliable components and cost is the cost function: 

Cy — Cj + Pj Cp 


where C T is the total cost, Cj is the manufacturing cost, P f is the failure probability, and C F is the cost 
incurred from structural failure. The manufacturing cost is represented by 
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( 9 ) 


Ci«£ Cjffp + q, 

j=l 

where N is the number of distribution parameters, Cj(Pj) is the manufacturing cost for the j* distri- 
bution parameter, and C 0 is constant cost. The total cost can be minimized when 

BCt 

-JL = o j=i, n ( 10 ) 

Substituting equation (10) into equation (9) achieves optimization by solving a system of nonlinear 
algebraic equations: 


3C j (p j> , r _ 34>(-P) Sp 

~Se~ F ~ 5jT" ^ 


0 


j=l, N 


(ID 


For a normally distributed random variable, 9p/9Pj can be calculated by equations (5) and (6). 


DEMONSTRATION FOR OPTIMIZATION OF SMART COMPOSITE STRUCTURES 
Maximization of Structural Reliability with Bounded Distribution Parameters 

Reliability maximization is demonstrated using an airfoil-like composite shell structure with an 
adaptive torque plate (fig. 3). The structure consists of skin, spar, and a torque plate. The skin and spar 
are a graphite-epoxy composite. The torque plate consists of five layers with the configuration 
(45/-45/steel/-45/45) in which the ±45° plies are control plies. In each control ply, both control 
(hybridizing actuation) and regular strips exist. However, in this paper, the control strip is assigned 
throughout the control ply for computational simplicity. The actuation material volume ratio is the per- 
centage of actuation composite material in a ply. The percentage of the actuation fiber in a composite 
is denoted by the actuation fiber volume ratio. The constituent material properties for regular plies and 
their assumed probabilistic distribution and coefficient of variation (representing the range of the 
scatter) are summarized in table I. The corresponding fabrication variables used to make the smart 
composite wing are summarized in table H, and those for the control are summarized in table in. 

The torque plate is fixed at the root. When control plies are electrically actuated, the angle of 
attack varies about its reference position. The assumed design criterion requires that the actuated 
change in the angle of attack be greater than 3°. The mean value m; of i* random variable is assumed 
to be bounded by 


0.97 m 0 < m ; < 1.03 nu 

u i 1 v i 


( 12 ) 


where niQ is i* reference mean. The standard deviation Oj of the i 111 random variable is assumed to be 
bounded by 
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0.6 o 0 < o t 


(13) 


where o 0 is the i* reference standard deviation. The actuated change in the angle of attack due to 
control is first probabilistically evaluated at the reference mean and standard deviation. The sensitivity 
factors for the six most sensitive random variables are presented in figure 4, which shows that the 
actuation material volume ratio has the largest (absolute) sensitivity factor (-0.568), followed by the 
actuation strain coefficient, the actuation material fiber modulus, the expansion and contraction electric 
field strengths, and the actuation material fiber volume ratio (-0.540, -0.409, -0.220, -0.217, -0.179, 
respectively). The reliability can be improved by increasing the mean value of these six random vari- 
ables as indicated by equation (5). Also, from equation (6), the smaller the standard deviation of the 
random variable, the larger the reliability of the structure. Therefore, the largest mean value and the 
smallest standard deviation should be used to achieve maximum reliability. For the purpose of demon- 
stration, only mean values and standard deviations of the actuation material volume ratio and actuation 
strain coefficient are chosen to be distribution parameters. The probability density function (pdf) of the 
actuated change in the angle of attack for a 3-percent increase in the mean of either the actuation 
material volume ratio or the actuation strain coefficient is shown in figure 5. Note that the mean value 
of the actuated angle increases whereas the scatter remains the same. The reliability can be improved 
by moving the pdf away from the design critical value. The pdf of the actuated angle for a 40-percent 
reduction in the scatter of either random variable is also shown in figure 5. The mean response (a 
change in the angle of attack) remains the same. The scatter of the response is reduced the most with a 
40-percent reduction in the scatter of the random variable with the largest sensitivity factor. The 
increase in the reliability index p for a 3-percent increase in the values of the actuation material 
volume ratio and the actuation strain coefficient are estimated using equation (5). The increase in p 
due to a 40-percent reduction in the scatter is estimated by equation (6). The estimated values of P 
shown in table IV are verified by reanalyzing the smart composite structure using IP ACS computer 
code. From the analysis, the maximum reliability index under the distribution parameter constraints is 
3.712, which corresponds to a 0.9999-reliability. 


Minimization of Reliability-Based Cost 

The cost minimization is demonstrated by using an airfoil-like composite shell structure with an 
adaptive skin (fig. 6). The structure consists of skin and spar. The spar is a graphite-epoxy composite. 
The skin is a ±45° adaptive intraply hybrid fiber composite for control. The upper and lower skins are 
electrically actuated to generate a change in the angle of attack. The reliability of the structure is cal- 
culated for an assumed performance criterion which requires that the actuated change in the angle of 
attack be greater than 1.2°. 

To simplify the optimization procedure, only the mean value of the actuation material volume ratio 
for the upper skin is chosen to be the design parameter. The initial cost C[(P]) is assumed to be 

Cj(P j) = PjC (14] 

where Pj is the mean value of the actuation material volume ratio, and C is a constant; C and C p 
are assumed to be 1 and 20, respectively. Substituting equation (14) into equation (11) achieves 
optimization by solving the following equation: 
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( 15 ) 


e aoep) ap_ = 0 
F ap 3Pj 

Equation (15) is used to find the minimum cost at a 0.91 -actuation material volume ratio. The 
result is verified by plotting the total cost and failure probability versus the mean actuation material 
volume ratio shown in figure 7. As can be seen, increasing the mean actuation material volume ratio 
results in a lower failure probability. The expected cost for failure PfC F is also reduced. At the opti- 
mum mean actuation material volume ratio (0.91), the manufacturing cost and expected cost due to 
failure are balanced. However, the reduction in the expected cost due to failure is less than the 
increase in the manufacturing cost when the mean actuation material volume ratio is greater than 0.91. 
This study demonstrates that a selection among possible arrangements for the optimum cost with a 
reliability consideration can be achieved using the probabilistic method adopted in the computer code 
IPACS. 


CONCLUSIONS 

Optimization problems with reliability considerations for adaptive intraply hybrid fiber composite 
structures were studied and demonstrated. The study was accomplished by integrating probabilistic 
methods and smart composite concepts into the IPACS (Integrated Probabilistic Assessment of Com- 
posite Structures) computer code. Probabilistic sensitivity factors from the probabilistic assessment of 
composite structures were used in the optimization procedure. The reliability for actuated change in the 
angle of attack of an airfoil-like composite shell structure with an adaptive torque plate was maximized 
to 0.9999 probability with constraints on the mean and standard deviation of the actuation material 
volume ratio (percentage of actuation composite material in a ply) and the actuation strain coefficient. 
This demonstrated that structural reliability with constraints on the distribution parameters could be 
maximized. The reliability-based cost for an airfoil-like composite shell structure with adaptive skin 
was also minimized with the mean actuation material volume ratio as the design parameter. The 
minimum cost was found when the mean actuation material volume ratio was equal to 0.91. Therefore, 
the balance between the maximum reliability and the minimum cost was achieved. 
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TABLE I.— STATISTICS OF FIBER AND 
MATRIX PROPERTIES FOR 
GRAPHITE-EPOXY COMPOSITE 


[Assumed distribution type, normal; assumed 
coefficient of variation 0.03.] 


Property 

Assumed 

mean 

Fiber modulus direction, Mpsi 


Longitudinal, E fll 

31.0 

Transverse, E^ 

2.0 

Fiber shear modulus, Mpsi 


IN-plane, G n2 

2.0 

Out-of-plane, G^j 

1.0 

Fiber Poisson’s ratio 


In-plane, v n2 

.2 

Out-of-plane, v 23 

.25 

Matrix 


Elastic modulus, E m , Mpsi 

.5 

Shear modulus, G m , Mpsi 

.185 

Poisson's ratio, v m 

.35 


TABLE n.— STATISTICS OF 
FABRICATION VARIABLES 
[Assumed distribution type, 
normal; assumed coeffi- 
cient of variation, 0.03.] 


Variable 

Assumed 

mean 

Volume ratio 


Fiber 

0.60 

Void 

.02 

Ply 


Misalignment angle, 3 deg 

0 

Thickness, in. 

.005 


a Assumed coefficient of variation, 
I (stdv). 
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TABLE m — STATISTICS OF CONTROL- 
RELATED VARIABLES 
[Assumed distribution type, normal; assumed 
coefficient of variation, 0.03.] 


Variable 

Mean 

Actuation material fiber modulus, 

12.4 

Mpsi 


Actuation fiber volume ratio 

0.60 

Actuation material volume ratio 

0.80 

Actuation strain coefficient, in TV 

10"* 

Electric field strength, V/in. 

10 5 

Expansion 

Contraction 

I0 5 


TABLE IV.— EHNHANCED RELIABILITY BY MANUFACTURING CONTROLLING RANDOM VARIABLES 1 


Random variable 
being controlled 
(actuation material) 

Sensitivity 

factor 

Reliability index, fi 

AP b 

P 

Ap c 

P 

Three-percent increase in mean 

Forty-percent reduction in scatter 

Equation (3) (estimation) 

IPACS 

Equation (4) (estimation) 

IPACS 

Actuation material 








Volume ratio 

-0.568 

0.568 

2.388 

2.349 

0.1 11 

1.258 

1.249 

Strain coefficient 

-.540 

.540 

2.360 

2.301 

.212 

2.032 

2.066 

Fiber volume ratio 

-0.179 







Fiber modulus 

-0.409 







Electric field strength 








Expansion 

-0.220 







Contraction 

-.217 







Initial (unchanged) 


— 

1.820 

1.820 




Final (combined) 


— 

3.375 

3.712 





“Estimated using sensitivity factors and recalculated by IP ACS. 

b Ap* = -iL Am-. 

onij 

■A|T = i!l 4o,. 
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Sensor/control strip expansion 
Sensor/control strip contraction 

Induced moment 


Figure 1 . — Adaptation of intraply hybrid to smart composite system, (a) Intraply hybrid composite system, (b) Structural control 
using sensor/control material. 
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Figure 2.— Concept of probabilistic assessment of composite structures. 
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Figure 3. — Airfoil-type composite shell structure with adaptive torque plate. All dimensions are in inches. 



Sensitivity factors 

Figure 4. — Sensitivity factors of uncertain variables to reliability for lower bound of change in angle of attack. 
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Probability density function 
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Actuation material volume ratio 
3-percent mean increase 
40-percent scatter reduction 
Actuation strain coefficient 
3-percent mean increase 
40-percent scatter reduction 


initial 

(unchanged) 

Final 

(combined) 


Red line 


Failure 



3 3.2 3.4 3.6 3.8 

Angle of attack, deg 

Figure 5. — Actuated change in angle of attack. 
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Figure 6. — Airfoil-type composite shell structure with adaptive skin. All dimensions are in inches. 
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Reliability-based total cost, Cj 



.84 .87 .90 .93 .96 .99 

Mean control volume ratio 

Figure 7. — Reliability-based total cost as function of mean control volume ratio. 
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